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ABSTRACT:  Theoretical  calculations  describing  the  Initiation 
in  heterogeneous  TNT,  RDX,  Tetryl,  Coop.  B,  Pentollte,  and 
75/25  CyclotoW'  produced  by  shocks  up  to  27  kllobars,  are 
given,  me  mot  spot*  Initiation  mechanism  Is  simulated  by 
using  appropriate  equations  of  state  of  the  explosives  and 
the  resultant  growth  from  shock  to  detonation  wave  Is  shown 
to  be  In  qualitative  agreement  with  experimental  results. 

The  Shock  wave  In  the  explosive  travels  with  Increasing 
velocity  due  to  the  release  of  energy  during  chemical  reaction 
In  the  neighborhood  of  the  shock  front,  and  then  takes  on  a 
constant  value  upon  reaching  full  detonation  velocity. 
Numerical  experiments  show  that  the  qualitative  ordering  of 
the  sensitivity  of  the  6  explosives,  based  on  one-dlmenslonal 
gap  test  calculations.  Is  the  same  as  that  obtained  from  the 
U.S,  Nival  Ordnance  Laboratory  experimental  gap  tests. 
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KQMR  62-160 
INTRODUCTION 

♦  The  processes  occurring  in  the  shock  initiation  of 

heterogeneous  solid  explosives  are  far  more  difficult 
to  understand  than  those  occurring  in  homogeneous  liquids. 
The  heterogeneous  nature  of  the  solid  leads  to  neighbor¬ 
hoods  of  small  local  convergences  in  the  shock  nave 
pattern  with  the  result  that  reaction  takes  place  in 
these  neighborhoods  first.  The  energies  released  at 
these  "hot  spots"  and  the  distribution  and  number  of 
these  hot  spots  then  serve  to  determine  the  initiation 
characteristics  of  the  explosive.  A  quantitative  theo- 
retlcal  calculation  taking  these  interactions  into  account 
,  presents  a  fluid  flow  problem  so  formidable  as  to  be 

beyond  the  abilities  of  present  methods  and  computers. 

Tet  by  lumping  these  interactions  together  so  that  the 
individual  actions  are  lost  but  the  qualitative  effect  of 
the  uhole  remains,  it  is  possible  to  derive  results  which 
are  in  excellent  qualitative  agreement  with  the  experimen¬ 
tal  results.  This  is  the  purpose  of  this  paper.  &  previous 
theoretical  paper1  has  discussed  the  initiation  of  liquid 
explosives  which  show  a  behavior  that  is  different  from 
that  of  solids. 

The  large  number  of  experimental  investigations  and 
4  more  detailed  explanations  of  the  associated  possible 

initiation  mechanisms  have  been  previously  reviewed  1  . 
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However,  two  experimental  techniques  will  bo  briefly 
described  hero  elnco  the  theoretical  results  reported 
in  this  paper  are  based  upon  a  Mathematical  Model  which 
is  taken  as  an  approximation  to  the  experiments. 

An  Important  explosive  sensitivity  test  Is  the 
Ngap  test1*.  Detailed  descriptions  of  the  test  as  used 
at  the  U.  S.  Xaval  Ordnance  Laboratory  are  given  else¬ 
where®'^1^  the  experimental  arrar^ement,  as  shown  In 
Pig.  1,  consists  of  a  cylindrical  donor  explosive  charge, 
a  gap  built  up  of  an  Inert  material,  the  cylindrical 
acceptor  explosive  charge  whose  sensitivity  Is  to  be 
determined,  a  small  air  gap,  and  a  test  plate.  The  donor 
(tetryl  pellets)  Is  caused  to  detonate  through  initiation 
by  a  detonator.  A  detonation  wave  propagates  through  the 
donor  and  subsequently  a  shock  wave  Is  transmitted  into 
the  gap.  This  shock  wave,  continuously  attenuated  by 
rarefaction  waves  from  behind  and  from  the  sides,  passes 
into  the  acceptor  explosive  and,  if  sufficiently  strong, 
creates  a  temperature  high  enough  to  cause  detonation. 

The  energy  stimulus  In  this  case  is  the  transmitted  shock. 
A  reaction  strength  comparable  to  detonation  is  indicated 
by  a  hole  punched  in  the  steel  test  plate.  The  donor 
charge  material  and  dimensions  are  kept  fixed  and  the 
length  of  the  gap  is  increased  or  decreased  until  that 
critical  length  d^Q,  called  the  " 50 %  gap  value”  or 
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"5Q0  point",  Is  reached  such  that  further  increase  or 
decrease  will  cause  the  acceptor  explosive  always  to 
fall  to  detonate  or  always  to  detonate,  respectively. 

These  5 00  gap  values,  under  carefully  controlled  condi¬ 
tions,  are  quite  reproducible  and  reaaxfcably  sharp; 
statistical  analysis  of  results  froa  a  50  shot  run  on 

deep  B  showed  that  the  500  point  did  not  vary  by  a ore 

N  11 

than  >  fraction  of  a  millimeter  .  On  this  sensitivity 

scale,  the  larger  the  500  gap  value  the  wore  sensitive 
the  explosive.  Saear  caaeras,  high  speed  framing  cameras, 
and  electrical  probes12,  and  most  recently,  a  conducting 

wire  along  the  axis  of  the  aaaoably13  hare  been  used  to 

\ 

study  the  build-up  froa  a  shock  wave  to  a  detonation  wave 

ll 

In  the  acceptor  explosive  .  The  pressure  of  the  initial 
transaltted  shock  In  the  acceptor  explosive,  that  corres¬ 
ponds  to  d^Q,  la  called  the  "initiating  pressure"  or  the 
"500  pressure",  P^q.  The  experimental  gap  test  Involves 
two  -dimensional  hydrodynamic  flow  since,  for  practical 
reasons,  the  characteristic  disaster  to  donor  +  gap  length 
ratio  Is  not  large  enough  to  rule  out  the  effect  of 
lateral  rarefaction  waves  on  the  amplitude  of  the  shock 
when  It  enters  the  acceptor  (even  along  the  axis  of 
symmetry). 

A  clever  experimental  technique,  which  minimizes 
troublesome  two-dimensional  effects  in  the  shock  inltia- 


4 


BOLTH  62-160 


♦ 


* 


A 

i 


% 


tlon  of  solid  explosives,  is  the  SQL  wedge  test  developed 
lA 

by  Jacobs  .  Plane  wave  lenses  are  oonblned  with  addi¬ 
tional  high  explosive  and  an  Inert  shock  attenuator  to 
produce  ~  chock  wave  of  desired  aaplitude  in  the  acceptor 
explosive.  Ytt  acceptor  explosive  in  the  shape  of  a  wedge 
is  Minted  as  shofcu  in  Pig.  2  so  that  the  progress  of  the 
shock  or  detonation  wave  can  be  seen  as  notion  along  the 
slant  face.  For  fixed  booster  geometry,  the  thickness  of 
the  Inert  attenuator  is  varied  until  it  is  possible  to 
discern  a  build-up  of  the  shock  velocity  in  the  acceptor 
to  detonation  velocity.  Various  workers  have  postulated 
that  this  build-tap  is  due  to  the  continuous  energy  release 
behind  the  shock  front  brought  about  by  chemical  reaction. 
This  technique  has  been  further  refined  and  used  in  the 
very  extensive  work  of  Cambell,  et  al1^.  Recent  use  of 
the  wedge  test  to  measure  initiating  pressures  for  Cong)  B-3 
has  resulted  in  values  somewhat  larger  (ca.  5096)  than 
those  obtained  in  the  gap  test  of  Fig.  l1^. 

The  mechanism  for  the  release  of  chemical  energy 
behind  the  shock  front  in  heterogeneous  explosives  is 
quite  complicated  and  still  poorly  understood.  But  it  is 
reasonable  to  attribute  the  difference  between  the  initia¬ 
tion  behaviors  of  homogeneous  and  heterogeneous  explosives 
to  the  presence  of  voids  and  other  defects  in  the  latter. 
Experimental  results  for  nltromethane*^  and  nitrone thane - 
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tlon  of  solid  explosives,  is  the  POL  wedge  test  developed 
i4 

by  Jacobs  .  Plane  ware  lenses  are  eonblned  with  addi¬ 
tional  high  explosive  and  an  inert  shock  attenuator  to 
produce  »  shock  wave  of  desired  amplitude  in  the  acceptor 
explosive.  kte  acceptor  explosive  in  the  shape  of  a  wedge 
is  mounted  as  shorn  in  Mg.  2  so  that  the  progress  of  the 
shock  or  detonation  wave  can  be  seen  as  notion  along  the 
slant  face.  For  fixed  booster  geometry ,  the  thickness  of 
the  Inert  attenuate r  is  varied  until  it  is  possible  to 
discern  a  build-up  of  the  shock  velocity  in  the  acceptor 
to  detonation  velocity.  Various  workers  have  postulated 
that  this  build-up  is  due  to  the  continuous  energy  release 
behind  the  shock  front  brought  about  by  chemical  reaction. 
This  technique  has  been  further  refined  and  used  In  the 
very  extensive  worlc  of  Campbell,  et  al1^.  Recent  use  of 
the  wedge  test  to  measure  initiating  pressures  for  Comp  B-3 
has  resulted  In  values  somewhat  larger  (ca.  20£)  than 
those  obtained  In  the  gap  test  of  Fig.  l1^. 

the  mechanism  for  the  release  of  chemical  energy 
behind  the  shock  front  in  heterogeneous  explosives  Is 
quite  complicated  and  still  poorly  understood.  But  It  Is 
reasonable  to  attribute  the  difference  between  the  Initia¬ 
tion  behaviors  of  homogeneous  and  heterogeneous  explosives 
to  the  presence  of  voids  and  other  defects  in  the  latter. 
Experimental  results  for  nltromethane1^  and  nitrone thane - 
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carborundum  mixtures1^  show  that  convergences  in  the 
mass  flow  and  impedance  mismatches  can  cause  local 
reactions  which  have  an  important  influence  on  the  initi¬ 
ation  process.  Per  the  mixture  a  detailed  knowledge  of 
the  shock  Interactions  is  more  important  than  a  knowledge 
of  the  values  of  the  thermochemical  constants.  It  is  well 
known  that  an  explosive  becomes  easier  to  initiate  as  the 
density  is  decreased;  and  it  is  concluded  that  the  Increase 
in  the  number  of  voids  leads  to  an  Increase  in  the  number  of 
hot  spots  formed,  either  through  reactions  occurring  on  the 

void  surface,  or  through  hot  spots  produced  by  shock  con- 

18 

vergence  beyond  the  void.  Recent  computer  calculations 
describe  the  fluid  flow  resulting  from  a  plane  shock  striking# 
from  below,  a  bubble  of  vacuum  suspended  in  nltromethane . 

The  bubble  is  closed  almost  simultaneously  everywhere 
because  the  lower  surface  has  almost  reversed  itself  by 
the  time  it  hits  the  upper.  The  collapse  is  accompanied  by 
the  generation  of  a  maximum  temperature  In  the  liquid,  just 
above  where  the  bubble  had  been;  the  temperature  generated 
is  at  least  twice  that  in  the  initial  shock. 

The  experimental  evidence  indicates  that  in  hetero¬ 
geneous  explosives  the  initial  shock,  which  is  too  weak 
to  raise  the  bulk  of  the  compressed  explosive  to  a  tempera¬ 
ture  sufficient  to  cause  reaction,  creates  hot  spots  where 
a  small  amount  of  energy  is  released  and  that  this  energy 
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production  ceases  soon  after  the  shock  front  has  passed 
over  the  particular  hot  spot.  Furthermore,  the  com¬ 
pression  waves  originating  at  these  centers  of  energy 
production  continuously  reinforce  the  shock,  leading  to 
eventual  detonation  at  or  near  the  front.  Among  others, 
Cachia  and  Whitbread1^  have  experimentally  determined  the 
conditions  under  which  the  initial  shock.  Instead  of 
accelerating  continuously  into  a  detonation  in  the 
acceptor,  decelerates  and  fades.  In  this  case  the  reac¬ 
tion  induced  by  the  shock  does  not  liberate  sufficient 
energy  to  overcome  the  losses  due  to  the  rarefaction  waves 
from  the  rear  and  sides.  Retonatlon  (detonation  backwards 
through  the  partially  reacted  explosive)  has  been  reported 
by  some  observers  but  not  by  others1^. 

In  this  paper  numerical  experiments,  based  on  a 
previously  discussed  computational  scheme1,  are  described 
in  which  a  hot  spot  mechanism  has  been  simulated  in  solid 
explosives  and  the  resultant  growth  from  shook  to  detona¬ 
tion  wave  is  shown  to  be  in  qualitative  agreement  with 
experimental  results.  The  hot  spots  are  simulated  by  an 
appropriate  choice  of  the  equation  of  state,  details  of 
which  are  given  in  Section  II.  In  Section  III,  the  criti¬ 
cal  gap  values  and  "50$  pressures"  derived  from  one- 
dimensional  gap  test  calculations  for  TNT,  RDX,  Tetryl, 
Comp  B,  Fentollte,  and  75/25  Cyclotol  are  compared  with 
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experimental  values  and  the  initiation  mechanism  i£? 
discussed  In  detail. 

II.  EQUATIONS  OP  STATE 

In  the  present  calculations,  the  equations  of  state 
of  unreacted  explosive  and  product  gas  are  similar  to 
those  previously  used1,  l.e.,  for  the  explosive, 

B,  -  *2  +  [(P+B)V,  -  (F°+B)Vj]  /(tb-1),  (1) 

*8  -  +  °r,»  1  +  [AVaI_T*  +  (y,-DW,/T,-(P0+B)v2]/(t,-1).  (2) 


*  -  K)T»  [p°+Vy,-(y,-D  «t,s*°''vS]  . 

(3) 

B  -  pSK)2, 

(3a) 

and  for  the  product  gas. 

V"  *v<vx,» 

W 

Kg  ”  av,eT’ 

(5) 

where  E,P,V,T,  cy,  c,  and  p*l/V>  are  respectively,  the 

specific  Internal  energy,  pressure,  specific  volume, 

teoperature,  specific  heat  at  constant  volume,  adiabatic 

sound  speed,  and  density;  and  r_  are  constants;  and  the 

e  6 

subscripts  g,  s,  and  0  refer  respectively  to  product  gas, 
unreacted  explosive,  and  the  initial  state.  If  Q  is  the 
change  In  specific  internal  energy  at  P°  and  T°,  l.e.. 
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0  -  *8(p0,i®)  -  *g(P°,*°),  (6) 

then  substitution  of  Bqs.  (2)  and  (5)  Into  (6)  defines 

Bja. 

K  -  *  +  °v,gT°*  (7) 

The  specific  energy  E  and  volume  V  for  a  mass  cell  of 
mixture  are  given  by 

E  -  wEg  +  (l-w)Eg,  V  -  Iffg  +  (l-w)Vg,  (8) 

where  w  (o  sc  w  *  l)  is  the  mass  fraction  of  unreacted 
explosive  in  the  cell.  The  solid  and  gas  in  the  cell  are 
assumed  to  have  the  same  instantaneous  values  of  P  and  T. 
All  the  constants  appearing  in  the  equations  of  state  for 
TUT,  RDX,  Tetryl,  Coop  B,  Pentolite,  and  75/23  Cyclotol 
are  listed  in  Table  I.  The  values  of  c°  for  the  6  explo¬ 
sives  were  obtained  from  experimental  shock  data  for  TNT  and 
Co*p  B16  by  interpolation  on  the  initial  densities  in 

or\  on 

the  Pfu  plane  v;  the  values  of  Y8  were  then  obtained 
from  experimental  values  of  P^q. 

While  the  same  form  of  the  equation  of  state  has 
been  used  for  the  solid  explosive  as  previously1  used  for 
the  liquid  TOT,  the  values  of  y  are  quite  different 
(y  •  3*178  for  liquid  TOT).  If  T^  is  the  temperature  of 
the  unreaoted  explosive,  then  the  temperature  along  the 
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lsentroplc  expansion  curve  from  the  state  T^  is  given 
*T 

*  -  T1(T1A) Vl  •  (9) 

Assuming  now  that  only  a  small  amount  of  taction  has 
already  occurred  before  expansion,  the  temperature  of  the 
mixture  can  still  be  approximated  by  £q.  (9).  For  solids, 
values  of  6.11  <  Yg  *  7*12  (see  Table  I)  lead  to  much 
greater  temperature  Increases  during  compression  and 
decreases  during  expansion  than  for  liquid  TNT  for  the 
name  relative  Increase  or  decrease  In  density.  The  result 
Is  that  a  rarefaction  wave  can  abruptly  quench  reaction  In 
a  cell  of  such  a  solid  explosive,  whereas  that  same  wave 
(mtme  in  the  sense  of  having  the  Identical  Instantaneous 
pressure-distance  curve  )may  not  quench  reaction  In  liquid 
TNT.  This  has  been  verified  In  numerical  experiments  for 
both  types  of  explosive. 

As  will  be  seen,  the  present  choice  of  yB  for  the 
solid  explosive  does  lead  to  effects  which  are  qualita¬ 
tively  Identical  with  those  produced  by  hot  spots.  This 
equation  of  state,  giving  bulk  temperatures  which  are  too 
high  for  the  compressions  encountered  In  the  gap  test, 
simulates  the  high  temperatures  at  the  hot  spots  (which 
are  surrounded  by  relatively  cold  explosive).  Thus,  the 
temperatures  calculated  here  may  be  considered  as  the 
"hot  spot  temperature s".  Likewise,  the  very  large  temperature 
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decrease  during  expansion  nay  be  considered  as  the  very 
large  decrease  in  "hot  spot  temperature"  due  to  the 
expansion  occurring  after  the  passage  of  the  shock. 

op 

The  inert  gap  material  is  assumed  to  have  a  P,V 
Hugoniot  curve  given  by 

P  -  58.4  (0.847-V6)/(0.272+Yo)2  , 

where  P  and  V„  are  In  kllobars  and  cor/g,  reepeetlvely; 
the  isentropic  expansion  curve  is  assumed  to  follow  the 
Hugoniot. 

III.  ONE-DIKENSIOMAL  GAP  TEST  CALCULATIONS 

The  calculations  described  below  assume  that  the 
fluid  flow  is  one-dimensional.  While  this  is  well  approxi¬ 
mated  by  the  experimental  wedge  test,  it  is  a  rougher  approx¬ 
imation  for  the  experimental  gap  test,  since  here  the 
pressure  amplitude  and  the  pressure  distribution  behind 
the  shock,  as  it  approaches  the  interface,  has  already  been 
Influenced  by  rarefaction  waves  from  the  sides. 

In  the  explosive  an  irreversible  first-order  chemical 
reaction. 

-  t«i  »  ud 

is  assumed,  where  the  chemical  kinetic  equation  which 
governs  the  conversion  of  unreacted  solid  explosive  [s]  to 
product  gas  [g]  is 

dw/dt  ■  -Z  w  exp  [-  E*  /  (RT)J  .  (12) 
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Here,  t,  E*,  Z,  and  R  are  respectively,  the  tine,  acti¬ 
vation  energy#  frequency  factor,  and  gas  constant.  Values 
for  E*  and  Z  are  given  in  Table  Z. 

The  calculatlonal  Method  involving  the  solution  of 
the  Lagrangean  hydrodynamic  equations,  the  equations  of 
state,  and  Eq.  (12)  by  finite  differences  is  described 
elsewhere1.  The  number  of  mesh  points  used  in  these 
computations  ranged  from  110  to  150  a 

The  one-dimensional  idealized  version  of  the  gap  test 
is  shown  in  Pig.  5.  A  number  of  numerical  experiments 
were  performed  in  which  the  length  of  the  donor  explosive 
(Tetryl),  dd<mor,  was  fixed  at  5.08  cm  and  the  length  of 
the  inert  gap  dg  was  varied  until  detonation  in  the  acceptor 
explosive  (TNT,  RDX,  Tetryl,  Comp  B,  Pentollte,  and  73/25 
Cydotol)  was  Just  barely  possible.  For  larger  gaps  than 
this  critical  length  d£,  detonation  never  propagated;  for 
smaller  gaps,  detonation  always  propagated.  The  results 
of  these  machine  computations  (on  an  IBM  7090  computer) 
are  given  in  Table  II.  Here  Pcrit  and  Tcrit  are  the 
values  of  the  pressure  and  temperature  of  the  transmitted 
shock  (as  it  entered  the  acceptor  explosive)  which  corres¬ 
pond  to  d£.  In  all  cases,  the  shock  initiated  some  reaction 
upon  entering  the  acceptor  explosive  and  had  advanced  a 
distance  xfl  into  the  explosive,  continuously  increasing  in 
amplitude2^,  at  the  time  the  explosive  mass  cell  just 
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FIG. 3  IDEALIZED  ONE- DIMENSIONAL  GAP  TEST 
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behind  the  shock  front  vent  to  complete  reaction,  the 
values  of  and  are  compared,  In  table  H,  with 

the  ejqperlmental  *50f  gap  values  d^  and  pressures  V 


obtained  at  the  TJ.S.  Maval  Ordnance  Laboratory  with  Luclte 
gaps11.  However,  the  Initial  densities  used  In  the  calcu¬ 
lations  and  those  used  In  the  experiments  are  different 
and  are  given  In  table  ZZ.  It  Is  encouraging  to  note 
that  although  Pcrlt  and  are  quite  different  for  some 
of  the  explosives,  the  relative  ordering  of  the  sensitivity 
of  the  explosives  on  the  P^^  scale  Is  the  sane  as  that 


on  the  F^q  scale 


at 


this  is  true  not  only  for  the  pure 
explosives,  MT,  RBX,  and  fetryl,  bnt  also  for  the  Mix¬ 
tures.  the  seam  ordering  holds  when  comparing  on  the  dj 
and  d^Q  scales. 

this  theoretical  nodel,  with  Its  numerous  approxima¬ 


tions,  can  be  expected  to  reproduce  experimental  trends 


for  series  of  pure  explosives  of  appreciable  differences 
in  sensitivity.’  the  fact  that  It  also  reproduces  trends 
for  the  mixtures  can  be  attributed,  in  part,  to  the  senrl- 
tlvity  of  the  results  to  the  choice  of  yg.  This  Is  shown 
most  clearly  for  RDX  and  Pentollte;  both  had  p°  *  1.59 


and  therefore,  by  the  aethod  used,  the  same  Bugoniot. 
However,  the  respective  ranges  of  0*P«r7  and  0*P«?11.8  kbar 
resulted  In  ym  values  of  7*12  and  6.99.  These,  in  turn, 

B 

resulted  In  the  same  Por^^  although  the  attenuation  (d^) 
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differed  bj  22  on  of  Inert  gap  mateiial.  Of  course,  am 

additional  difficulty  in  this  case  «u  the  greater  uncer- 

* 

talaty  of  the  proper  raises  for  E  ,  Z,  end  Q  of  Pentolite. 

The  same  sensitivity  to  y  Is  evident  in  the  series  of 
Hf4  SS3K  and  their  Mixtures,  75/25  Cyclotol  and  Coqp  B. 
Certainly  It  Is  to  he  expected  that  the  cyclotol  sitfa  75 % 
of  the  more  sensitive  component  sould  he  more  sensitive 
than  Coop  S  (about  Sof  HEX).  Both  experhaental  and  ooepu- 
tatlonal  results  reverse  the  expected  order*  Tbs  reason 
for  the  reversal  of  the  experimental  vaults  Is  probably 
lack  of  adequate  control  of  physical  variables  such  as 
particle  alee  of  the  IS,  Indeed,  recant  carefully  con¬ 
trolled  tests  of  pressed  charges  have  exhibited  the  expected 
sensitivity  ordering1^.  The  reversal  of  the  computational 
results  probably  arises  from  the  ranges  used  to  derive  y 
and  these  In  turn  are  based  05  the  experimental  results. 

It  seems  quite  likely  that  repeating  the  treatment  with  use 
of  the  experimental  results  for  pressed  charges  mould 
reverse  the  order  obtained  with  the  cast  changes  reported 
here. 

There  seems  little  doubt  that  the  present  treatment  can 

reproduce  experimental  trends  for  TNT,  Tetryl,  RDX,  and 

their  mixtures;  a  plot  of  P QT ^  vs.  Is  a  smooth  curve. 

The  departure  cf  the  Pentollte  result  from  this  curve  Is 

* 

attributed  chiefly  to  the  less  adequate  values  of  E  and  Z 
for  PSTN. 
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The  pressure  distributions  at  various  times  for  75/25 
Cyclotol  are  shown  In  Fig*  4;  the  other  explosives  show 
essentially  the  same  profiles.  In  Fig.  4(a) ,  the  Taylor 
detonation  wave  has  already  reflected  off  the  product 
gas  -  Inert  gap  Interface  and  a  reflected  rarefaction 
wave  Is  moving  hack  Into  the  product  gas,  and  a  transmitted 
shock  is  moving  Into  the  gap.  Only  the  rightmost  part  of 
the  gas  Is  shown  since  most  of  it  is  expanding  to  the  left 
Into  a  vacuum  (approximating  air  at  several  bars).  In 
Fig.  4(h)  the  progress  of  the  shook  in  the  gap  Is  shown 
and  In  Fig.  4(c)  the  shock  has  reflected  off  the  gap  - 
75/25  Cyclotol  Interface  with  the  result  that  higher  ampli¬ 
tude  shocks  are  moving  hack  Into  the  gap  and  Into  the 
explosive.  In  Figs.  4(d)  and  4(e),  It  can  he  seen  that 
the  pressure  (and,  therefor,  temperature)  of  the  trans¬ 
mitted  shock  Is  Increasing  because  of  chemical  reaction 
behind  the  front  (the  reaction  is  apparent  In  the  values 
of  w<l  In  the  numerical  solutions)  although  the  rarefaction 
wave  from  the  rear  would  tend  to  attenuate  the  shock.  The 
pressure  (and  temperature)  near  the  Interface  continuously 
decreases  and,  with  this  decrease,  the  chemical  reaction 
abruptly  ceases  In  contrast  to  the  liquid  TMT  case1  (where 
the  temperature  continues  to  rise).  Here  the  temperature 
rises  monotonically  from  the  interface  to  the  shock  front 
and  only  in  the  neighborhood  of  the  front  is  there 
appreciable  chemical  reaction.  In  Figs.  4(f)  and  4(g),  the 
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Fig.  4  The  pressure  distribution  in  the  one- 

dimensional  gap  test  for  various  values 
of  t  for  75/25  Cyclotol.  The  heavy 
vertical  lines  A,  B,  C,  and  D  represent 
respectively,  the  air-product  gas  inter¬ 
face,  the  attenuator-75/25  Cyclotol  (or 
product  gas)  interface,  and  the  75/25 
Cyclotol- air  Interface.  Each  dot  repre¬ 
sents  a  mass  point  whose  initial  position 
is  given  by  the  ijgrangean  coordinate  x. 

The  motion  of  any  mass  point  in  the  X 
direction  can  be  followed  by  observing 
the  motion  of  that  same  dot  (provided 
the  observer  has  a  pair  of  telescopic  eyes). 
The  values  of  t  (tisec),  as  measured  from 
the  instant  at  which  the  shock  first 
entered  the  explosive,  are  as  follows: 

(»)  -108.4,  (b)  -25.3,  (0)  -4.1,  (d)  0.4, 

(•)  4.8,  (f)  9.1.  («)  13.5,  (h)  16.0, 

(1)  17.5,  (J)  25.3. 


20 


IBISSIBBSBi 

HanssHBy 

piiiSiij 

nSpanranpl 
■■■■■■■!■ M 

■■■■■■■■■ll 

HIHIIIIIll 

Nniinuil 

iiiiiliily 

SSSSSSSKS] 

■■■■■■■■■SI 


NOLTR  62-160 


reaction  is  sufficiently  rapid  so  as  to  lead  to  a  rapid 
increase  in  pressure  near  the  shock  front.  The  transition 
to  detonation  Is  completed  by  Fig.  4(h).  The  detonation  and 
retonatlon  waves  are  shown  In  Fig.  4(1).  Retonatlon  com¬ 
pleted,  is  shown  in  Fig.  *+(j). 

The  t,X  diagram  is  particularly  Informative  for 
showing  the  wave  and  particle  motion  and  this  is  depicted 
In  Fig.  5  for  75/25  Cyclotol.  Only  motion  in  the  vicinity 
of  the  gap-explosive  Interface  Is  shown.  The  curve 
represents  the  motion  of  the  shock  wave  in  the  gap  before 
It  strikes  the  gap-explosive  Interface.  Cg  represents 
the  path  of  the  interface  which  moves  with  the  local  par¬ 
ticle  velocity  while  a  transmitted  shock  wave  moves  along 
Into  the  explosive.  Increasing  In  velocity  from  A  to  B. 
High-order  detonation  occurs  just  behind  the  shock  front 
at  B,  and  a  constant  velocity  detonation  wave  advances 
along  Cjp  A  retonatlon  wave2^  starting  at  B  moves  through 
the  partially  reacted  pre-compressed  explosive  along 
until  it  contacts  the  Interface  at  D,  where  It  gives  rise 
to  a  transmitted  shock  that  moves  Into  the  gap  along  Cg. 
Since  the  Instantaneous  temperature  decreases  monotonlcally 
from  the  shook  front  to  the  Interface  In  the  partially 
reacted  pre-compressed  explosive  (and  the  temperature  at 
the  Interface  also  decreases  with  time  after  passage  of 
the  shock).  It  is  possible  to  find  a  curve  which  Is  the 

locus  of  points  at  which  reaction  is  barely  significant. 
Therefore,  along  lines  of  constant  t  to  the  left  of  Cj  the 
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reaction  baa  essentially  ceased  (because  T  is  too  low), 
and  to  the  right  of  Cj  reaction  is  significant.  The 
lengths  of  th.  horizontal  linos  In  the  region  between  Cj 
and  Cj  represent  the  instantaneous  widths  of  the  reaction 
zones  which  supply  an  Increasing  amount  of  energy  to  drive 
the  shock  front  along  C y  Thus,  the  model  used  for  the 
above  computations  gives  results  which  are  in  qualitative 
agreement  with  the  hot  spot  theory  of  shock  initiation  of 
solid  explosives.  These  hot  spots  are  considered  to  be 
"active"  and  feeding  chemical  energy  to  the  shock  wave 
only  for  a  short  time  after  the  shook  wave  has  passed  over 
them  (and  before  the  rarefaction  wave  from  the  rear  has 
sufficiently  cooled  them).  The  hot  spots  may  be  considered 
to  be  distributed  in  the  reaction  zone  that  lies  between 
and  C^.  These  time-dependent  reaction  zones  are  not 
to  be  confused  with  the  steady  state  reaction  zones;  the 
latter,  of  course,  are  smaller. 
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Let  the  subscripts  1  and  2  refer  TKT  and  Comp  B,  respec¬ 
tively.  Prom  Ref.  16 

UjL  *  c^  +  a^u,  i  *  1,2  (PI) 

on  their  respective  Hugoniots,  where  n  and  u  are  the 
shock  and  particle  velocities  respectively,  and  the  c^ 
and  are  constants.  For  the  initial  densities 
p^«  1.615  g/cm^  and  p^  -  1.72  g/cm^,  c^  ■  0.259 
cm/p-sec,  c^  •  0.271  cm/|isec,  ■  2.05,  and  a2  **  1.86 
Consider  a  test  explosive  of  initial  density  p°  with  an 
equation  of  state  given  by  Bq.  (1),  and  assume  that  at 
fixed  u  the  shock  pressure  P  can  be  found  by  a  linear 
interpolation  (or  extrapolation)  on  the  initial  densities 
of  TNT  and  Coop  B,  i.e., 

P  -  Pj  -  (P1-P2)(p(l)-p°)/(p(l)-p(S)).  (F2) 

where 

Px  -  p(1)(c{1)+a1u)u,  1  -  1,2  (F3) 

is  obtained  from  Eq.  (PI)  and  the  Hugoniot  relation, 

P  -  pM  Uu,  J  -  0,  1,  2  .  (F4) 
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the  Initial  pressure  F°  is  always  neglected  in  these 
calculations  and  the  subscript  s  has  been  dropped. 

Substitution  of  Bq.  (F2)  in  (F2)  leads  to 

P  -  p°(c°+  an)u,  (P5) 

Where 

p°c°  -[pWcW(pW  "  pc,V1)c(1)(p«2>.p0)]/(p(1).  p^)t 
P°«  -jp^2^«2(p^1^-p°)-p^1^»i(p^2^-p0)j/(p^1^-p^s^)*  (*6) 

Kllnlnatlng  P  betwam  Eqa.  (M)  and(F5)  gives 

D  -  o°  +  au  (P7) 

for  the  test  explosive.  Thus  cj  -  o°  is  taovn  in  Eq.  (ja). 
Values  of  c®  and  a  are  given  in  Table  I. 

9 

Eliminating  E,  P,  and  V  between  Eqs.  (1),  (Ft)  for  J-O,  and 
the  other  two  Hugoniot  relations, 

B  -  E°  -  i  P(V°-V),  (P8) 

u2  -  P(V°-V>,  (P9) 

leads  to 

T  -  1  -  2  [l  -  ^  +  (e°)2/(Uu)J  ,  (FiO) 

Where  the  subscript  s  has  been  dropped.  Obviously  Bq. 

(FIO)  does  not  relate  V  linearly  to  u  as  demanded  by 
Bq.  (F7)  if  y  -  constant.  Therefore  y  must  be  found  such 
that  tne  U,u  curves  of  Bqs.  (77)  and  (FIO)  are  reasonably 
dose.  Substituting  Eq.  (F7)  into  (FIO)  gives 

y  -  1  +  2ja-l  +  ac°/(c°+au)J,  * 

l.e.,  y  depends  on  u.  From  Eq,  (F5)  it  is  found  that 


(Fll) 
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2  *b  -  c°(l+4*P/B)^  ,  (F12) 

where  the  positive  root  has  been  taken.  Using  this  value 
in  Xq.  (m)  jri.lda 

riP)  -  142  {»-l4a«/C  (l-MtaP/B)*  4  13}  (F13) 

The  values  of  a  and  B  given  in  Table  Z  and  the  experimental 
m3Cf  pressures"  listed  In  Table  ZZ  are  used  to  evaluate 
an  average  y,  y_w  ,  over  the  Bugoniot  pressure  Interval, 

•T  • 

0«P<dP^Q,  by  means  of  the  relation, 

ya=  r„.  -  i  [t(o)^(p50)]  .  (*i4) 

Bor  the  six  solid  explosives  oonsldered,  y(0)  and  t(p^q)  vary 
by  at  most  120  (TNT  value)  from  yB,  which  is  considered 

suri  iciently  accurate  considering  the  approximations  which 
have  been  made  already. 

22  See  footnote  14  of  Ref.  (l) . 

2^  A  similar  Increase  In  shock  save  amplitude  was  obtained  by 
N.  H.  Boyer  (and  others),  Aeronutronlc  Systems  Inc.  (Ford 
Motor  Co.),  Pub.  Ho.  U-500,  November  15,  1958,  who  used  the 

O  y 

very  low  values,  E*  «  46  cal/g  and  Z  *  10  sec’*,  which  were 

calculated  from  the  equation  representing  the  adiabatic 

explosion  time  of  a  homogeneous  explosive  and  the  experimental 

14 

delay  times  of  Majowlcs  and  Jacobs  for  heterogeneous  solid 
explosives.  These  low  values  greatly  reduced  the  temperature 
sensitivity  of  the  reaction  rate. 


2$ 
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^  ?crlt  caaa  oa3l^3r  nade  to  coincide  wit*  P^q  by 
appropriately  adjusting  the  value  of  Ta,  cy,  £*,  or  Z, 
but  then  the  purpose  of  the  computations  would  be  lost* 

25  Caspbell,  et  al^  find  no  retonation  wave  in  their  expert, 
meats  and  explain  this  by  postulating  a  change  in  the  hot 
spot  characteristics  of  the  partially  reacted  pre- 
compressed  explosive.  The  model  used  in  the  numerical 
calculations  cannot  take  these  changes  into  consideration 
and  therefore  yields  a  retonation  wave. 
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